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Abstract 

^^ I Pade Approximants can be used to go beyond Vector Meson Dominance in 

a systematic approximation. We illustrate this fact with the case of the pion 
vector form factor and extract values for the first two coefficients of its Taylor 
expansion. Pade Approximants are shown to be a useful and simple tool for 

^ ' incorporating high-energy information, allowing an improved determination of 

^ . these Taylor coefficients. 
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1 Introduction 

It has been known for a long time that the pion vector form factor (VFF) in the space- 
hke region is very well described by a monopole ansatz of the type given by Vector 
Meson Dominance (VMD) in terms of the rho meson. However, it has remained unclear 
whether there is a good reason for this from QCD or it is just a mere coincidence and, 
consequently, it is not known how to go about improving on this ansatz. 

To begin our discussion, let us define the form factor, F((5^), by the matrix element 

(7r+(p')| I u^u - ^ d^d - \ sYs |7r+(p)) = {p + pY F{Q') , (1) 

where Q^ = —{p' — p)'^, such that Q^ > corresponds to space-like data. Since the 
spectral function for the corresponding dispersive integral for F^Q^) starts at twice the 
pion mass, the form factor can be approximated by a Taylor expansion in powers of 
the momentum for \Q'^\ < (2m7r)^. At low momentum, Chiral Perturbation Theory is 
the best tool for organizing the pion interaction in a systematic expansion in powers of 
momenta and quark masses [H [21 E] ■ With every order in the expansion, there comes 
a new set of coupling constants, the so-called low-energy constants (LECs), which 
encode all the QCD physics from higher energies. This means, in particular, that the 
coefficients in the Taylor expansion can be expressed in terms of these LECs and powers 
of the quark masses. Consequently, by learning about the low-energy expansion, one 
may indirectly extract important information about QCD. 

In principle, the coefficients in the Taylor expansion may be obtained by means of 
a polynomial fit to the experimental data in the space-like region |j below Q^ = 4m^. 
However, such a polynomial fit implies a tradeoff. Although, in order to decrease 
the (systematic) error of the truncated Taylor expansion, it is clearly better to go 
to a low-momentum region, this also downsizes the set of data points included in 
the fit which, in turn, increases the (statistical) error. In order to achieve a smaller 
statistical error one would have to include experimental data from higher energies, i.e. 
from Q^ > Ami. Since this is not possible in a polynomial fit, the use of alternative 
mathematical descriptions may be a better strategy. 

One such description, which includes time-like data as well, is based on the use 
of the Roy equations and Omnes dispersion relations. This is the avenue followed 
by [H E], which has already produced interesting results on the scalar channel ^, 
and which can also be applied to the vector channel. Other procedures have relied 
on conformal transformations for the joint analysis of both time-like and space-like 
data [7j, or subtracted Omnes relations [SIE]- Further analyses may be found in Ref. 

m- 

On the other hand, as already mentioned above, one may also consider an ansatz 
of the type 

m\... = (i + -j^) ' • (2) 

V VMD / 



^Time-like data is provided by tttt production experiments and, consequently, they necessarily 
correspond to values of the momentum above the tttt cut, i.e. \Q^\ > 4toJ with Q^ < 0. 



Even though the simphcity of the form of Eq. ([2]) is quite astonishing, it reproduces 
the space-hke data rather well, even for a range of momentum of the order of a few 
GeV, i.e. Q^ ^ 4m^. If this fact is not merely a fluke, it could certainly be interesting 
to consider the form ([2]) as the first step in a systematic approximation, which would 
then allow improvement on this VMD ansatz. 

In this article, we would like to point out that the previous VMD ansatz for the 
form factor ^ can be viewed as the first element in a sequence of Pade Approximants 
(PAs) which can be constructed in a systematic way. By considering higher-order terms 
in the sequence, one may be able to describe the space-hke data with an increasing 
level of accuracy q Of course, whether this is actually the case and the sequence is 
a convergent one in the strict mathematical sense or, on the contrary, the sequence 
eventually diverges, remains to be seen. But the important difference with respect 
to the traditional VMD approach is that, as a Pade sequence, the approximation is 
well-defined and can be systematically improved upon. 

Although polynomial fitting is more common, in general, rational approximants 
(i.e. ratios of two polynomials) are able to approximate the original function in a 
much broader range in momentum than a polynomial [TT]. This will be the great 
advantage of the Pades compared to other methods: they allow the inclusion of low 
and high energy information in a rather simple way which, furthermore, can in principle 
be systematically improved upon. In certain cases, like when the form factor obeys a 
dispersion relation given in terms of a positive definite spectral function (i.e. becomes a 
Stieltjes function), it is known that the Pade sequence is convergent everywhere on the 
complex plane, except on the physical cut [T2]. Another case of particular interest is in 
the limit of an infinite number of colors in which the form factor becomes a meromorphic 
function. In this case there is also a theorem which guarantees convergence of the Pade 
sequence everywhere in a compact region of the complex plane, except perhaps at a 
finite number of points (which include the poles in the spectrum contained in that 
region) [13J. In the real world, in which a general form factor has a complicated 
analytic structure with a cut, and whose spectral function is not positive definite, we 
do not know of any mathematical result assuring the convergence of a Pade sequence 
|14j . One just has to try the approximation on the data to learn what happens. 

In this work we have found that, to the precision allowed by the experimental 
data, there are sequences of PAs which improve on the lowest order VMD result in a 
rather systematic way. This has allowed us to extract the values of the lowest-order 
coefficients of the low-energy expansion. 

We would like to emphasize that, strictly speaking, the Pade Approximants to a 
given function are ratios of two polynomials Pn{z) and Qm{z) (with degree N and M , 
respectively), constructed such that the Taylor expansion around the origin exactly 
coincides with that of f{z) up to the highest possible order, i.e. f{z) — R^i^) = 
0{z^'^~^'^~^^). However, in our case the Taylor coefficients are not known. They are, in 
fact, the information we are seeking for. Our strategy will consist in determining these 
coefficients by a least-squares fit of a Pade Approximant to the vector form factor data 



^Obviously, unlike the space-like data, one should not expect to reproduce the time-like data since 
a Pade Approximant contains only isolated poles and cannot reproduce a time-like cut. 



in the space-like region. 

There are several types of PAs that may be considered. In order to achieve a 
fast numerical convergence, the choice of which one to use is largely determined by 
the analytic properties of the function to be approximated. In this regard, a glance 
at the time-like data of the pion form factor makes it obvious that the form factor 
is clearly dominated by the rho meson contribution. The effect of higher resonance 
states, although present, is much more suppressed. In these circumstances the natural 
choice is a P^ Pade sequence [TT], i.e. the ratio of a polynomial of degree L over a 
polynomial of degree oncQ Notice that, from this perspective, the VMD ansatz in (j2j) 
is nothing but the P° Pade Approximant. 

However, to test the aforementioned single-pole dominance, one should check the 
degree to which the contribution from resonances other than the rho may be neglected. 
Consequently, we have also considered the sequence P2, and the results confirm those 
found with the PAs P^. Furthermore, for completeness, we have also considered the so- 
called Pade- Type approximants (PTs) [151 116]. These are rational approximants whose 
poles are predetermined at some fixed values, which we take to be the physical masses 
since they are known. Notice that this is different from the case of the ordinary PAs, 
whose poles are left free and determined by the fit. Finally, we have also considered 
an intermediate case, the so-called Partial- Pade approximants (PPs) [15], in which 
some of the poles are predetermined (again at the physical masses) and some are left 
free. We have fitted all these versions of rational approximants to all the available 
pion VFF space-like data [l7]-[22]. The result of the fit is rather independent of the 
kind of rational approximant sequence used and all the results show consistency among 
themselves. 

The structure of this letter is the following. In section [2] we begin by testing the 
efficiency of the P/" Fades with the help of a model. In section[3]we apply this very same 
method to the experimental VFF. Firstly, in sec. 13. H we use the Fade Approximants 
P/"; then, in Sec. 13. 2^ this result is cross-checked with a P^ PA. Finally, in sec. 13. 3^ we 
study the Fade- Type and Partial-Fade approximants. The outcome of these analyses 
is combined in section H] and some conclusions are extracted. 

2 A warm-up model 

In order to illustrate the usefulness of the PAs as fitting functions in the way we propose 
here, we will first use a phenomenological model as a theoretical laboratory to check 
our method. Furthermore, the model will also give us an idea about the size of possible 
systematic uncertainties. 

We will consider a VFF phase-shift with the right threshold behavior and with 
approximately the physical values of the rho mass and width. The form-factor is 



■'Conventionally, without loss of generality, the polynomial in the denominator is normalized to 
unity at the origin. 



recovered through a once-subtracted Omnes relation, 



F{Q' 



exp 



Q2 roo 



dt 



m 



tit + g2 



(3) 



where S{t), which plays the role of the vector form factor phase- shift [HI [9l [23], is given 
by 

^Mptpit) 



6{t) = tan"^ 
with the t-dependent width given by 



Mj-t 



(4) 






e{t- 4ml) 



(5) 



and a{t) = a/1 — Arh'^/t. The input parameters are chosen to be close to their physical 
values: 



0.15 GeV 



M. 



0.6 GeV' 



4m: = 0.1 GeV' 



(6) 



We emphasize that the model defined by the expressions (JSJIS]) should be considered 
as quite realistic. In fact, it has been used in Ref. [H [9l [23] for extracting the values 
for the physical mass and width of the rho meson through a direct fit to the (timelike) 
experimental data. 

Expanding F{Q^) in Eq. ([3]) in powers of Q^ we readily obtain 



F{Q' 



1 - ai Q2 + a2 Q' - as Q' + 



(7) 



with known values for the coefficients a^. In what follows, we will use Eq. (jTj) as the 
definition of the coefficients Oj. To try to recreate the situation of the experimental 
data p!7]-[22] with the model, we have generated fifty "data" points in the region 
0.01 < Q^ < 0.25, thirty data points in the interval 0.25 < Q' < 3, and seven points 
for 3 < Q' < 10 (all these momenta in units of GeV'). These points are taken with 
vanishing error bars since our purpose here is to estimate the systematic error derived 
purely from our approximate description of the form factor. 

We have fitted a sequence of Pade Approximants P/'(Q') to these data points and, 
upon expansion of the Fades around Q' = 0, we have used them to predict the values 
of the coefficients Oj. The comparison may be found in Table [H The last PA we have 
fitted to these data is Pf . Notice that the pole position of the Fades differs from the 
true mass of the model, given in Eq. ([6]). 

A quick look at Table [1] shows that the sequence seems to converge to the exact 
result, although in a hierarchical way, i.e. much faster for ai than for 02, and this one 
much faster than 03, etc... The relative error achieved in determining the coefficients a^ 
by the last Fade, Pf , is respectively 0.9%, 8% and 26% for ai, 02 and 03. Naively, one 
would expect these results to improve as the resonance width decreases since the P^ 
contains only a simple pole, and this is indeed what happens. Repeating this exercise 
with the model, but with a Fq = 0.015 GeV (10 times smaller than the previous one). 





A" 


PI 


Pf 


Pf 


Pt 


Pf 


Pf 


F(Q2) (exact) 


ai(GeV-2) 
a2 (GeV-4) 
a3(GeV-6) 


1.549 
2.399 
3.717 


1.615 
2.679 

4.444 


1.639 
2.809 
4.823 


1.651 
2.892 
5.097 


1.660 
2.967 
5.368 


1.665 
3.020 

5.579 


1.670 
3.074 

5.817 


1.685 
3.331 

7.898 


Sp(GeV2) 


0.646 


0.603 


0.582 


0.567 


0.552 


0.540 


0.526 


0.6 



Table 1: Results of the various fits to the form factor F(Q'^) in the model, Eq. ([3|). The 
exact values for the coefficients Oj in Eq. ([7]) are given on the last column. The last row 
shows the predictions for the corresponding pole for each Fade (sp), to be compared to the 
true mass M^ =0.6 GeV^ in the model. 

the relative error achieved by Pf for the same coefficients as before is 0.12%, 1.1% 
and 4.7%. On the other hand, a model with Tq five times bigger than the first one 
produces, respectively, differences of 2.1%, 14.4% and 37.8%. 

As we have mentioned in the introduction, it is possible to build a variation of the 
PAs, the Pade-Type Approximants, where one fixes the pole in the denominator at the 
physical mass and only the numerator is fitted. We have also studied the convergence 
of this kind of rational approximant with the model. Thus, in this case, we have 
placed the P^ pole at Sp = M^ and found a similar pattern as in Tabled] For Pf, the 
Pade-Type coefficient ai differs a 2.5% from its exact value, 02 by 16% and 03 by 40%. 

Based on the previous results, we will take the values in Table [Has a rough estimate 
of the systematic uncertainties when fitting to the experimental data in the following 
sections. Since, as we will see, the best fit to the experimental data comes from the 
Pade P/, we will take the error in Tabled] from this Pade as a reasonable estimate and 
add to the final error an extra systematic uncertainty of 1.5% and 10% for ai and 02 
(respectively). 

3 The pion vector form factor 

We will use all the available experimental data in the space-like region, which may 
be found in Refs. [l7]-[22]. These data range in momentum from Q^ = 0.015 up to 
10 GeV^. 

As discussed in the introduction, the prominent role of the rho meson contribution 
motivates that we start with the Pi Pade sequence. 



3.1 Pade Approximants P^ 

Without any loss of generahty, a P/' Pade is given by 

L-l 



pHQ' 



1 + ^afe(-g2)fc + (_g 



2\L 



k=0 



ol 






l + ^^TT^Q"^ 



(8) 



where the vector current conservation condition Pi{0) = 1 has been imposed and the 
coefficients a^ are the low-energy coefficients of the corresponding Taylor expansion of 



Q^F(Q^) (GeVO 
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Q\Ge\) 



Figure 1: The sequence of Pf' PAs is compared to the available space-like data [I7j-[22]: P° 
(brown dashed), Pi (green thick-dashed), P^ (orange dot-dashed), Pf (blue long-dashed), 
Pf (red solid). 



the VFF (compare with (I7l) for the case of the model in the previous section). 

The fit of Pi to the space-like data points in Refs. [l7j-[22] determines the coeffi- 
cients ttfc that best interpolate them. According to Ref . [21] , the form factor is supposed 
to fall off as l/Q"^ (up to logarithms) at large values of Q^. This means that, for any 
value of L, one may expect to obtain a good fit only up to a finite value of Q^, but not 
for asymptotically large momentum. This is clearly seen in Fig. [H where the Fade 
sequence P(" is compared to the data up to L = 4. 

Fig. [2] shows the evolution of the fit results for the Taylor coefficients ai and 02 for 
the Pi PA from L = up to L = 4. As one can see, after a few Fades these coefficients 
become stable. Since the experimental data have non zero error it is only possible to 
fit a Pi PA up to a certain value for L. From this order on, the large error bars in 
the highest coefficient in the numerator polynomial make it compatible with zero and, 
therefore, it no longer makes sense to talk about a new element in the sequence. For 
the data in Refs. [IZ]-[22], this happened at L = 4 and this is why our plots stop 
at this value. Therefore, from the PA Pi we obtain our best fit and, upon expansion 
around Q^ = 0, this yields 



ai = 1.92 ± 0.03 GeV^ 



aa = 3.49 ± 0.26 GeV" 



(9) 



with a xVdof = 117/90. 

Eq. ([8]) shows that the pole of each Pi PA is determined by the ratio Sp = ai/ai+i- 
This ratio is shown in Fig. [31 together with a gray band whose width is given by iMpFp 
for comparison. From this figure one can see that the position of the pole of the PA is 
close to the physical value M^ [25], although it does not necessarily agree with it, as 
we already saw in the model of the previous section. 
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Figure 2: ai and 02 Taylor coefficients for the Pf' PA sequence. 
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Figure 3: Position Sp of the pole for the different P^. The range with the physical values 
Mp ± MpTp is shown (gray band) for comparison. 

3.2 Cominent on P2 Pades 

Although the tiine-Hke data of the pion form factor is clearly dominated by the p(770) 
contribution, consideration of two-pole P2 PAs will give us a way to assess any possible 
systematic bias in our previous analysis, which was limited to only single-pole PAs. 

We have found that the results of the fits of P2 PAs to the data tend to reproduce 
the VMD pattern found for the P/" PAs in the previous section. The P2 PAs place the 
first of the two poles around the rho mass, while the second wanders around the complex 
momentum plane together with a close-by zero in the numerator. This association of a 
pole and a close-by zero is what is called a "defect" in the mathematical literature |26j. 
A defect is only a local perturbation and, at any finite distance from it, its effect is 
essentially negligible. This has the net effect that the P^ Pade in the euclidean region 
looks just like a P^ approximant and, therefore, yields essentially the same results. 
For example, for the P2, one gets 



ai = 1.924 ± 0.029 GeV~ 



a2 = 3.50 ±0.14 GeV" 



(10) 



with a xVdof = 120/92. 



3.3 Pade Type and Partial Pade Approximants 

Besides the ordinary Pade Approximants one may consider other kinds of rational 
approximants. These are the Pade Type and Partial Pade Approximants [131 [13, [IS] • 
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Figure 4: Low-energy coefficient ai from the T^ Fade- Type sequence 

In the Fade Type Approximants (FTAs) the poles of the Fade are fixed to certain 
particular values, which in this context are naturally the physical masses. On the 
other hand, in the Fartial Fade Approximants (FFAs) one has an intermediate situation 
beytween the FAs and the FTAs in which some poles are fixed while others are left as 
free parameters to ffi. 

Since the value of the physical rho mass is known {Mp = 775.5 MeV), it is natural 
to attempt a fit of FTAs to the data with a pole fixed at that mass. The corresponding 
sequence will be called T^. This has the obvious advantage that the number of param- 
eters in the fit decreases by one and allows one to go a little further in the sequence. 
Our best value is then given by the Fade Type Approximant Tf, whose expansion 
around Q^ = yields the following values for the Taylor coefficients: 

ai = 1.90 ± 0.03 GeV~2 , aa = 3.28 ± 0.09 GeV"^ , (11) 

with a xVdof = 118/90. 

The previous analysis of FTAs may be extended by making further use of our 
knowledge of the vector spectroscopy |^2_5i|. For instance, by taking Mp = 775.5 MeV, 
Mpi = 1459 MeV and Mpn = 1720 MeVp we may construct further Fade- Type se- 
quences of the form T2 and Tg^. 

In the FTA sequence T2 one needs to provide the value of two poles. For the ffist 
pole, the natural choice is M^. For the second pole, we found that choosing either M^ 
or M^„ (the second vector excitation) does not make any difference. Both outcomes 
are compared in Fig. ([5]) . Using M^, , we found that the T2 FTA yields the best values 
as 

ai = 1.902 ± 0.024 GeV-^ as = 3.29 ± 0.07 GeV-\ (12) 

with a xVdof = 118/92. 

Using M^„ as the second pole one also gets the best value from the T| FTA, with 
the following results: 

ai = 1.899 ± 0.023 GeV"2^ as = 3.27 ± 0.06 GeV-^ (13) 



*As will be seen, results do not depend on the precise value chosen for these masses. 
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Figure 5: Low energy coefficient ai for tlie Tg" Pade-Type sequence witli Mp and Mpi (left), 
and with Mp and Mpn (riglit). 

with a x^/dof = 119/92. We find the stability of the results for the coefficients ai 2 
quite reassuring. 

We have also performed an analysis of the PTA sequence T^, with similar conclu- 
sions. From the T| we obtain the following values for the coefficients: 



ai = 1.904 ±0.023 GeV" 



a2 = 3.29 ± 0.09 GeV" 



(14) 



with a xVdof = 119/92. 

Finally, to complete our analysis, we will also consider Partial Pade Approximants, 
in which only part of the denominator is given in advance. In particular, we study the 
PPA sequence P^^ |j in which the first pole is given by M^ and the other is left free. 
The best determination of the Taylor coefficients is given by Pii, and they yield 

ai = 1.902 ±0.029 GeV"^^ 02 = 3.28 ± 0.09 GeV-^ (15) 

with the free pole of the PPA given by M]^^^ = (1.6±0.4 GeV)^ and a xV^of = 119/92. 

4 Combined Results and conclusions 

Combining all the previous rational approximants results in an average given by 
a, = 1.907 ±0.010stat±0.03syst GeV"' , 02 = 3.30 ± 0.03stat ± 0.33syst GeV-^ (16) 

The first error comes from combining the results of the different fits by means of a 
weighted average. On top of that, we have added what we believe to be a conservative 
estimate of the theoretical (i.e. systematic) error based on the analysis of the VFF 
model in Sec. [21 We expect the latter to give an estimate for the systematic uncertainty 
due to the approximation of the physical form factor with rational functions. For 
comparison with previous analyses, we also provide in Table [2]the value of the quadratic 
radius, which is given by (r^) = 6 ai . 



^See Ref. [13] for notation. 





ir') (fm^) 


02 (GeV-^) 


This work 


0.445 ± 0.002,tat ± O.OOTsyst 


3.30 ± 0.03stat ± 0.33syst 


CGL [31^5] 


0.435 ±0.005 




TY ^ 


0.432 ±0.001 


3.84 ±0.02 


BCT [27] 


0.437 ±0.016 


3.85 ±0.60 


PP [9J 


0.430 ±0.012 


3.79 ±0.04 


Lattice 28 


0.418 ±0.031 





Table 2: Our results for the quadratic radius (r^) and second derivative 02 are compared 
to other determinations [H El [3 [271 El HH]. Our first error is statistical. The second one is 
systematic, based on the analysis of the VFF model in section 2. 



In summary, in this work we have used rational approximants as a tool for fitting 
the pion vector form factor. Because these approximants are capable of describing the 
region of large momentum, they may be better suited than polynomials for a description 
of the space-like data. As our results in Table 2 show, the errors achieved with these 
approximants are competitive with previous analyses existing in the literature. 
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